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The non-classical rotational inertia fraction of the identical cylindrical solid ^He below 300 mK is 
studied at 496 and 1173 Hz by a double resonance torsional oscillator. Below 35 mK, the fractions 
are the same at sufficiently low rim velocities. Above 35 mK, the fraction is greater for the higher 
than the lower mode. The dissipation peak of the lower mode occurs at a temperature ~ 4 mK lower 
than that of the higher mode. The drive dependence of the two modes shows that the reduction of 
the fraction is characterized by critical velocity, not amplitude nor acceleration. 
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Ever since the observation of "Non-Classical Rotation 
Inertia" (NCRI) of solid ''He [l| contained in their tor- 
sional oscillators by Kim and Chan 0,0]) apparent super- 
fluidity in solid "^He has motivated numerous experimen- 
tal and theoretical studies. A clear understanding of the 
observed apparent partial decoupling between solid ^He 
and the container walls of torsional oscillators is lack- 
ing at present. A number of interpretations of the ob- 
served decoupling has been introduced including glassy 
state [4], vortex liquid state [1], superglass Q, and a 
phenomenological two component model (see [1] for a 
recent review). Many of these theoretical ideas contain 
predictions on frequency dependence of the magnitude of 
decoupling. Although the frequency of previous torsional 
oscillators has been varied between 180 and 1500 Hz, the 
oscillators contained different independently grown solid 
^He samples in different cells and laboratories. Here, 
we report on our measurements of NCRI of the identi- 
cal solid "^He sample contained in a torsional oscillator 
having two resonant frequencies differing by a factor of 
2.4. We observe that the frequency dependence is a func- 
tion of both temperature (T) and the imposed oscillation 
velocity. However, at sufRciently low temperature and os- 
cillation velocity, the measured superfluid fraction does 
not depend on frequency. At T > 35 mK, the superfluid 
fraction depends on frequency, oscillation velocity and 
temperature. In addition to these frequency dependent 
effects, we show clear evidence for hysteretic behavior 
of NCRI phenomenon depending on the history of oscil- 
lation velocity and thermal processing from the normal 
state above 300 mK to low temperature. Our observa- 
tions put severe constraints on some of the theoretical 
notions on the interpretation of NCRI. 

The conventional torsional oscillator having one res- 
onant frequency is made of one solid ''He container at- 
tached to a single torsion rod. Our torsional oscillator as 
shown in Fig. [T]contains two masses, the upper "dummy" 
without solid *He and the lower containing sample ^He, 
and two torsion rods. The torsion rods, a portion of the 
dummy mass and the lower flange fitted to the sample 
container are all machined from a single stock of BeCu 
rod. The first resonant mode (frequency /i = 495.8 Hz) 
occurs when the two masses move in phase and the second 
mode (/2 = 1172.8 Hz) when the two masses move 180° 
out of phase. The phase relation of the modes are verified 



by measuring the oscillator response at both the upper 
and the lower detectors. Measurements of the oscillator 
response are made with the upper electrode fins. The 
lower electrode is used to calibrate oscillation amplitude 
at the sample cell against the upper detector response. 
Owing to the uncertainty in the capacitance between the 
lower fixed and movable electrodes, the absolute oscilla- 
tor amplitude is accurate only to within ±10 %. The 
relative accuracy in comparing the oscillator amplitudes 
between the two modes is ±5 %. The temperature at the 
mixing chamber is measured by a calibrated '^He melting 
pressure thermometer. The temperature of the copper 
isolator block which supports the oscillator is monitored 
by a carbon resistance thermometer. 




FIG. 1: Double resonance torsional oscillator. It contains 
two masses and BeCu torsion rods (a) (2 mm outer diame- 
ter). The upper mass includes a BeCu disc (b) with movable 
detector and generator electrode fins (c). The lower mass 
is made of a Stycast 1266 epoxy block (e) containing a hol- 
low cylindrical solid *He sample space (10.2 mm in diameter 
and 7.6 mm in height. A second movable detector electrode 
(f) attached to the bottom side of the lower mass is used to 
calibrate the upper electrode detector motion. *He is filled 
through a central 0.8 mm diameter hole (d) drilled in the tor- 
sion rods and the upper mass. Fixed electrodes placed facing 
the movable electrodes are not shown. 

The "background" characteristics [Q] of the resonant 
frequencies, fio (i = 1, 2), and the quality factors Qio of 
our torsional oscillator with the sample cell empty were 
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measured as functions of temperature and drive level. 
The Q,o (1-4 x 10^ and 0.63 x 10^ for i = 1 and 2, respec- 
tively, at 30 mK) was determined from the ring down 
time. The oscillation amplitude was calibrated against 
the measured Qio by the ring down time. The tempera- 
ture dependence of the resonant frequency and the qual- 
ity factor of the empty cell were both slightly dependent 
on the drive level but qualitatively similar to those ob- 
served by others. After the sample cell is loaded with 
solid ^He at 37 bar, for example, the resonant frequen- 
cies measured at 300 mK decreased from the empty cell 
values by A/i = 0.667 Hz and A/2 = 2.049 Hz for the 
first and second mode, respectively. 

Our sample solid ^He is grown from commercially avail- 
able helium gas with nominal '^He concentration of 0.5 
ppm by the blocked capillary method. The cell is loaded 
near 3 K to a desired pressure. The part of the fill line 
which is thermally anchored to the 1 K pot is quickly 
plugged with solid while ''He in the cell attached to the 
mixing chamber remains in the liquid state. The liquid 
pressures and temperatures, when the plug is formed and 
when the ''He cools down to the liquid/solid coexistence 
point, are monitored by a pressure sensor at the mixing 
chamber. The change in pressure within the cell itself is 
sensed from the oscillation amplitude of one of the modes 
during solidification. After the plug is formed, the total 
time elapsed to freeze ^He in the cell was varied from 
several minutes to four hours. 

When our torsional oscillator with the sample cell filled 
with solid ^He is cooled, the two resonant frequencies 
[fiiT)) increase below 300 mK and reach maxima near 
35 mK. This increase in frequency has been interpreted 
to arise from a partial decoupling of solid "^He from the 
sample cell wall and has been called NCRI [11, • The 
NCRI fraction (NCRIf) is evaluated from our measure- 
ments as [/.(T) ~ (/,;o(r) - A/,)] /A/,. The tempera- 
ture dependent changes in resonant frequencies are ac- 
companied by decrease in quality factor. A measure of 
dissipation induced by the presence of solid "^He is evalu- 
ated as 5[l/Qi{T)] = l/Q^{T) - l/Qm{T). The data of 
resonant frequency and quality factor reported here are 
obtained after stabilizing the temperature and waiting 
for transient effects to diminish to negligible levels. 

In all, we have grown and studied five solid samples 
with final solid pressures between 27 and 42 bar. The 
observed features were qualitatively similar in all solid 
samples. Among the samples we studied, we present re- 
sults, as an example, of a 37 bar sample solidified over 
four hours. The temperature of this sample was main- 
tained below 310 mK throughout all of the measurements 
reported here after the initial cool down so as to avoid any 
changes in the properties of our sample by some anneal- 
ing process. The measured NCRIf and the dissipation 
while keeping the rim velocity sufficiently low to avoid 
critical velocity effects (see below) are shown in Fig. [2] 
It should first be noted that the maximum NCRIf ob- 
served is only 0.1 %, which is considerably smaller than 
those found generally in annular cells [1, [llj . (The max- 
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FIG. 2: Temperature dependence of non-classical rotational 
inertia fraction and change of dissipation (see text) at /i (Q) 
and /2 (A). Pressure of the solid is 37 bar. Rim velocities 
are 8.5 ~ 13.5 /im/sec at /i, and 7.2 ~ 16.0 /^m/sec at /2. 



imum NCRIf in another sample solidified deliberately 
rapidly over several minutes was larger at 0.12 %.) On 
the other hand, the maximum is not so different from 
some of those found in cylindrical cells in other labora- 
tories [ll,[ll,[ll. At T < 35 mK, NCRIf is independent 
of both temperature and frequency. The dissipation is 
small in the same temperature range. In the higher tem- 
perature range T > 35 mK, both the NCRIf and the dis- 
sipation show significant dependence on frequency. The 
NCRIf measured with the first mode is lower than that 
with the second mode where T > 35 mK. At T = 100 mK, 
NCRIf(/2)/NCRIf(/i) = 1.9 and this ratio continues to 
increase at higher temperatures. The observed magni- 
tude of 6Q^^ is also smaller than in other annular cells. 
In the temperature range, 35 mK < T < 100 mK, SQ^^ 
is greater than <5Q^^ contrary to what might be expected 
from the relative difference of NCRIf measured with the 
two modes. The dissipation maximum value of the first 
mode is greater than that of the second by a factor of 
1.7. The temperature where the dissipation maximum 
occurs is slightly (^ 4 mK) lower in the first mode than 
the second. The temperature "width" of the dissipation 
peak is narrower in the first mode. 

All previous torsional oscillator experiments have 
shown that the NCRIf decreases as the drive level is in- 
creased and the decrease has been ascribed as a "critical 
velocity" effect. It has not been clarified if the reduc- 
tion in NCRIf was controlled instead by the amplitude of 
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displacement or acceleration. To investigate the critical 
velocity effect, the oscillator is cooled down to 19 mK in 
the first mode with a low velocity near 10 /im/s. When 
the drive level is increased, to our surprise, the NCRIf 
does not diminish but remains constant up to 600 ^m/s 
as shown by circles in Fig. [3l It is conceivable that the 
NCRIf decays if we wait long enough, but the time con- 
stant for this is estimated (assuming exponential decay 
for this purpose) to be more than 10^ hours. Following 
the advice by A. Clark and M. Chan [l3|, the torsional 
oscillation is initiated in the first mode with a relatively 
high drive level at 300 mK. While keeping the drive level 
constant, the oscillator is then cooled down to 19 mK, 
where the rim velocity becomes 610 /im/s (shown by 
crosses). In this process, as the drive level is decreased, 
the NCRIf increases and eventually attains the low veloc- 
ity limit consistent with Fig. [2l This increase in NCRIf 
as the rim velocity is decreased is similar to those velocity 
effects as observed by others [^,^3, 13, 15]. If the oscilla- 
tion is initiated at any other rim velocity < 610 /xm/s at 
300 mK and cooled down to 19 mK, the observed NCRIf 
is that shown by the cross in Fig. [3] at the correspond- 
ing rim velocity. When similar procedures are followed 
with the second mode, the same NCRIf is obtained when 
plotted as a function of rim velocity (shown by trian- 
gles and diamonds). Cooling the oscillator at a low or 
high rim velocity is analogous to "zero field cooling" or 
"field cooling" process, as suggested to us by A. Clark 
and M. Chan [l2|, in the studies of superconducting and 
magnetic materials. 

The reduction in NCRIf as the drive level is increased 
is observed to be the same for both modes when it is 
plotted against the cell rim velocity. If NCRIf is plot- 
ted against displacement amplitude or acceleration, the 
reduction does not coincide in the two modes. This obser- 
vation demonstrates convincingly that it is the rim veloc- 
ity, not displacement amplitude (strain) nor acceleration 
(stress) applied to the solid, that best parameterize the 
reduction in NCRIf. The observed history dependence 
of NCRIf on the initial state set by the rim velocity at 
low temperature is an important characteristics of the 
supersolid state. 

Our measurements of critical velocity effects at 63 mK 
are shown in Fig. U) The circles (crosses) and diamonds 
(triangles) show NCRIf when the rim velocity is increased 
(decreased) for the first and second mode, respectively. 
The NCRIf decreases as the rim velocity is increased. 
When the difference at this temperature at low velocity 
(see Fig. [2]) is added to the NCRIf(/i), the reduction is 
best parameterized by the rim velocity just as at 19 mK. 
The saturation of the NCRIf below 15 fim/s as seen in 
Fig. [3] is absent. There is no history dependence at 63 
mK! The measured NCRIf does not depend on how the 
initial state is reached. The border between the history 
dependent low temperature behavior and the reversible 
higher temperature behavior appears to be around 40 
mK. This is close to the temperature where the NCRIf 
begins to decrease. 
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FIG. 3: Rim velocity dependence of non-classical rotational 
inertia fraction for two frequencies at the temperature of 19 
mK. Measurements are taken as the drive is decreased (-I-) 
(increased (O)) in sequence from the maximum (minimum) 
value shown for /i. Similarly, measurements are taken as the 
drive is decreased (A) (increased (o)) in sequence from the 
maximum (minimum) value shown for /2. Arrows indicate 
the direction of the rim velocity changes. 
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FIG. 4: Velocity dependence of non-classical rotational inertia 
fraction for two frequencies at 63 mK. Symbols refer to the 
same processes as in Fig. O 



Anderson interprets the observed NCRI as a conse- 
quence of vortex liquid and supercurrents flowing around 
a thermally excited fluctuation of vortices [1]. He sug- 
gests that there is a relaxation rate at which vortices can 
move in and out of the sample. The rate is thought to 
decrease as the temperature is lowered and as the ■^He 
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impurity concentration is increased. ^He is thought to 
act promote pinning. When the torsional osciUator fre- 
quency matches the rate of vortex motion, a maximum in 
dissipation would be observed. According to the model, 
the dissipation peak should occur at higher temperature 
for higher oscillator frequency. Our data shows that the 
dissipation peak for the second mode occurs at a slightly 
higher (by 4 ± 1 mK) temperature than that of the first 
mode. The direction of change in the peak temperature 
is consistent with the Anderson model, but a quantita- 
tive comparison is not possible (see Fig. 1 of Ref. 
If the Anderson model is correct, the vortex flow rate is 
greater than 500 s"'^ at T > 50 mK. It would then be 
expected that the NCRI is reversible, as indeed seen in 
our experiment, at the time scale (of order 10^ s) of our 
measurements. At low temperatures, T < 35 mK, how- 
ever, the observed relaxation rate is dependent on the 
history of velocity variation as illustrated in Fig. [3l The 
vortex motion rate is apparently dependent strongly on 
the existing number of vortices. 

From the point of view of vortex liquid model, the ir- 
reversible increase in NCRIf as seen in Fig. [3] from the 
high to low rim velocity may be interpreted as a process 
in which vortices can easily escape presumably from the 
cell surface. Once the vortices are absent, or nearly ab- 
sent, they are apparently excluded from entering solid 
^He even when the rim velocity is increased. These ob- 
servations are reminiscent of Meissner effect in supercon- 
ductors. A theoretical discussion of vortex nucleation in 
solid ^He has been given by Saslow ^Ifi] . 

According to the phenomenological two component 
model developed by Huse and Khandker [tI, if a solid 
sample is homogeneous, the change (Sf) in the torsional 
oscillator frequency in the T = limit is related to 
the maximum damping at the apparent supersolid tran- 
sition as Sf/f ~ S{l/Q)max- Evaluating the quanti- 
ties from the data shown in Fig. [21 we find the ratio 
{Sf / f)/S{l/Q)max to be 2.7 and 6.5 for the first and sec- 



ond mode, respectively. The ratios are comparable to 
those of the "best" among the samples studied by Kim 
and Chan jlB]- One might expect that the ratio should 
be independent of frequency if it is simply a measure of 
homogeneity. The temperature width of the dissipation 
peak is also thought to indicate the degree of sample 
inhomogeneity. It is apparent in Fig. [2] that the second 
mode shows a broader temperature width in 5(1/ Q) than 
the first. 

Nussinov et al. Q proposed that the observed increase 
in frequency and the peak in dissipation of torsional oscil- 
lators containing solid ^He might be explained in terms of 
solidification process of liquid-like component into a glass 
at low temperature. In their model of the glass state, 
there is a characteristic relaxation time s which increases 
as temperature is decreased. When the inverse relaxation 
time matches the oscillator frequency, a peak in the os- 
cillator dissipation appears. The change in frequency- 
dependent dissipation in their model naturally leads to 
a monotonic increase in frequency as temperature is de- 
creased. If s is assumed to take the simple form [3], 
s = soe^/*^^-^, where so and A are parameters, the mea- 
sured dissipations and frequency shifts of the two modes 
cannot be reconciled. Other forms of s might reduce dis- 
crepancy. Our observations place severe constraints on 
the glass model to explain the origin of observed NCRI 
effects of solid ^He. 
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